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ABSTRACT 



A study was made to relate certain optical properties 
with other observed properties of water sampled in southern 
Monterey Bay, California. Dominant wavelength, percent 
purity, and visual efficiency were determined for 6 5 near- 
shore water samples using a one-meter sample cell in a 
modified Beckman DU-2 spectrophotometer. Measurements made 
at the sample locations included salinity, surface temper- 
ature, phosphate, coliform count, oxygen, and particle size 
distribution, in the 1.04y to 27 .6y diameter range. Most 
of the sampling was done at or near the Monterey sewage 
outfall. Dominate wavelengths were found to vary between 
520 nm and 585 nm. Percent purity was found to fluctuate 
between 2 and 40 percent. Neither variable seems to be 
strongly sensitive to variation in treated sewage concen- 
tration for the Monterey outfall. The maximum transmission 
was linearly fitted to the logarithm of the projected 
particle area. A, by the equation: log-^(A/10 ) = .544 - .634 
(max T) , where A is in y , and transmission is per meter. 
Dominant wavelength for each of the 21 Forel-Ule scale colors 
was measured spectrophotometrically and compared with the 
dominant wavelengths of the samples. 
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I . INTRODUCTION 



A. BACKGROUND 

The effect of pollution on coastal waters is being 
increasingly studied. Of special interest are the dilution, 
the distribution, and the biological effects of effluents 
from sewage outfalls. Indications of sewage effluent dis- 
tribution and concentration are helpful in investigating 
the effect an outfall has on the marine ecology or the 
suitability for human recreational use of an area. Such 
studies require that quantitative measurements be made of 
the near-shore water. 

For many years scientists have observed the "color" and 

* 

the "transparency" of the world's oceans and lakes. These 
apparent properties vary with the different observers and 
conditions of observation [1] . Observations normally have 
not been made close to shore. Near-shore water is often in 
certain respects obviously different from other water in 
that its color tends to be "browner" and it contains more 
particulate matter. This is especially true in regions 
where effluents are being added from the shore. 



% 

"Color" or Forel-Ule color is obtained by comparing a 
series of numbered vials of colored solutions with water 
provided with a background (usually a Secchi disk lowered 
to a depth of one meter) . "Transparency" or Secchi depth 
is the depth (in meters) at which a Secchi disk (usually 
30 cm in diameter) disappears from view when lowered into 
the water. 
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Jerlov [2, pp . 145-150] has reported oceanic color on a 
more quantitative basis by using values calculated from 
spectral irradiance and radiance data. The emphasis, 
though, was on optically characterizing world oceanic 
waters . 

This study was pursued in an effort to obtain quantita- 
tive relationships between color signatures obtained from 
spectral transmittance (an inherent property) and other 
observed sea water properties. It was further hoped to 
check the practicality of determining effluent distributions 
from sewage outfalls by in situ (as opposed to remote) 
optical and electronic particle counting methods. 

B. NATURE OF THE PROBLEM 

1 . Optical Properties of Water 

In the study of the oceans one property of the water 
which has often been measured is its "color." The word 
"color" is an elusive term, since it depends upon the visual 
acuity of the observer and many other factors, such as sky- 
lighting, sea state, direction of observation, height of 
observer above the sea, and the nature of the reflecting 
background lowered into the water. In 1895 F.A. Forel put 
water color on a semi -quantitative basis with the introduc- 
tion of his color scale (31. Most of his measurements were 
made in Swiss lakes where the water was quite blue. The 
scale was extended by Ule, among others, to make possible 
the characterization of the more brownish waters found in 
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near-shore regions . A standard Forel-Ule scale is 
presently used for ocean water color determination [4, 

pp. 13 - 151 . 

A white surface (usually a Secchi disk) of 30-cm 
diameter is lowered into the ocean to a depth of one meter. 
This furnishes a reflecting background for the water, the 
color of which is matched with that of one of a series of 
vials containing colored solutions. The closest matching 
color yields the Forel-Ule number of the water. Addition- 
ally, the "transparency" of the water may be obtained. A 
Secchi disk is lowered in the water until it disappears and 
is raised until it is again sighted. The average of these 
two depths in meters is known as the "Secchi depth" or 
"transparency" [4, pp. 13-14]. 

The difficulty with both the Secchi depth and the Forel- 
Ule scale value is that they are at best semi-quantitative 
and depend rather strongly upon the visual acuity of the 
observer plus many other uncontrolled factors such as those 
suggested above. 

The measurement of inherent optical properties of sea 
water should allow samples from different depths and loca- 
tions to be more readily compared. One such property is 
that of beam transmittance , T, defined as the ratio of the 
intensity, I , of a beam of light of a particular wavelength 
having traveled a distance, d, in a water sample, to the 
intensity, I Q , of the light entering the sample. Beam 
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transmittance is related to the total beam attenuation 



coefficient , c, in the following way: 




Beam transmittance in water is highly dependent upon 
the wavelength of the incident energy. For pure water the 
maximum transmittance in the visible region occurs at about 
462 nm. As the water purity decreases, i.e. as the number 
of particles and energy absorbing impurities increase, the 
beam transmittance decreases and the point of maximum trans- 
mittance may shift. Dominant wavelength (color) and percent 
purity (purity of that color) of transmission for a given 
water sample can be assigned to that sample if the transmit- 
tance is known as a function of wavelength. 

The theory for the assignment of dominant wavelength 
(color) to samples is based upon the theory of trichro- 
maticity, the mechanism by which the human eye determines 
color. When the eye observes an image, the color is deter- 
mined by the mind's mixing of three independent colors, such 
as brown, orange, purple, etc. Thus, the mind has expressed 
a color or dominant wavelength which has three components: 
red, green and blue [2, pp . 140-142] . 

This theory can be put on a more quantitative basis by 
physically measuring the amounts of red, green,, and blue in 
an object or sample. The amount of energy in each of these 
three components is referred to as X, Y, and Z, respectively, 
the tristimulus values . Chromaticity coordinates 



10 



(trichromatic coefficients) x, y, and z, are simply the 
normalized tristimulus values. The energy required for each 
tristimulus value in order to produce a certain color has 
been calculated for various energy sources, such as sunlight 
(Source C) or a gas-filled tungsten lamp of color tempera- 
ture 2,850°K (Source A) (5, pp. 25^-316]. From the tristim- 
ulus values, chromaticity coordinates have been calculated. 
Figure 1 shows such a plot of chromaticity coordinates and 
the wavelength (color) each represents for Source A. The 
color-mixture envelope is shown as the outermost envelope 
on the graph [5, p. 2931 • For Source A (A, Fig. 1^ the 
chromaticity diagram has the distribution shown. If a dif- 
ferent source were used, the source point on the graDh would 
shift, and all lines of dominant wavelength would originate 
at the new "center" and extend to the color envelope. 

From measurments of the transmittance of light at various 
wavelengths in the visible region, the chromaticity coordi- 
nates, x and y, for transmission may be calculated for a 
sample. Once the coordinates are plotted on a chromaticity 
diagram, a line is drawn from the coordinates of the light 
source. A, through the sample point. The intersection of 
the color-mixture envelope and this line, W, marks the 
dominant wavelength of the sample. For point P in Figure 1, 
this is the point W. At point W on the color-mixture line 
the value of the wavelength is 595 nm. For point P the 



All figures are given in numerical order beginning on 
page 35 below. 



dominant wavelength is thus 595 nm. The percent purity of 
point P is defined as the ratio of the length of line seg- 
ment CP to that of segment CW, multiplied by 100. The value 
of the tristimulus value, Y, when multiplied by 100, is 
called the percent visual efficiency , a measure of the 
brightness or intensity of the energy passing through the 
sample . 

Once defined, dominant wavelength, percent purity, and 
visual efficiency of transmission are inherent properties 
of water; they completely express the color, purity of that 
color, and intensity of light transmission in the water. 
Additionally, this calculation reduces the values of trans- 
mittances observed throughout the visible spectrum to three 
numbers . 

2 . Near-shore Effluent Distributions 

For many years a standard measure of the effect of 
sewage outfall has been the coliform level. Coliform bac- 
teria are present in the digestive tracts of warm-blooded 
animals. By mapping the relative bacterial density, a kind 
of measure of the limits of a sewage distribution in the 
ocean may be obtained [6] . 

This indicator has become less reliable since many 
sewage treatment facilities now chlorinate their treated 
sewage effluent. The chlorine destroys most bacteria, 
including coliform, and makes the latter unreliable as an 
indicator of the distribution and dilution of sewage 
effluents . 
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The development of new and more accurate Indicators, or 
signatures, would aid greatly in marine ecological investi- 
gations, outfall design and location, and in the determina- 
tion of the suitability of coastal water for water-contact 
recreational use. 

C. METHOD OF INVESTIGATION 

A better indication of the distribution of sewage in 
the ocean may be obtained from the level of a particular 
nutrient, since even treated sewage is nutrient rich. 

An obvious difference between "normal" near-shore water 
and that near a sewage outfall may be its apparent or 
inherent color and the fact that there is often much addi- 
tional particulate matter present. On a quantitative basis 
these two parameters can provide indications of effluent 
distribution . 

The water of southern Monterey Bay provides a suitable 
region for the study of sewage effluent distributions. The 
bay is relatively sheltered and a number of outfalls extend 
from shore. Due to its accessibility, the Monterey sewage 
outfall has been the source of numerous studies both at the 
Naval Postgraduate School and at Hopkins Marine Station. 

The Monterey outfall extends 260 m into the bay. At its 
mouth the water depth is approximately 6 meters. The 
average flow rate from the outfall is 2.8 mgd of secondarily 
treated sewage effluent. 

To obtain an idea as to the measurable quantities which 
might be good effluent indicators, four stations were 
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selected around the Monterey sewage outfall for sampling. 
Each station was sampled at the surface and one meter from 
the bottom. Surface samples were also collected on one 
occasion near the Pacific Grove and Seaside sewage outfalls. 

To put optical data on a quantitative basis, dominant 
wavelength, percent purity, and visual efficiency were 
determined for each water sample. The data to compute these 
values were measured with a modified Beckman DU-2 spectro- 
photometer. Most laboratory spectrophotometers are equipped 
to handle cells from 1 to 10 cm in length. For such cells, 
measured transmittances for water are very nearly unity. 

Work with natural water requires the modification of such 
laboratory instruments to allow the light beam to pass 
through a relatively long sample. 

Particle distributions (size and number) were determined 
with an electronic particle counter. In addition, water 
temperature and samples for phosphate, salinity, oxygen, 
and coliform determination were taken. 

The data were then examined both graphically and statis- 
tically to determine the relationships, if any, between 
observables. The dominant wavelength of each of the 21 
standard Forel-Ule scale colors was also measured, as that 
information was not available. Such determinations make it 
possible to assign quantitative wavelength values to each 
of the solutions of the Forel-Ule scale. The dominant 






Model T Coulter counter having 15 channels. 



wavelength of the samples can then be related (indirectly) 
to a Forel-Ule color value. 
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II. OBSERVATIONAL PROGRAM 



A. SAMPLING LOCATIONS 

1 . Station Locations 

Samples were collected on six different days: 
September 16, 1971; September 28, 1971; October 5> 1971; 
October 12, 1971; October 18, 1971; and November 1 6 , 

1971. On October 18 stations -2 and -3 were sampled 

at the surface 350 m from their respective outfall. For 
detailed station locations see Figures 2 and 3 and 
Table I. 

a. Station M-l 

A surface sample was collected at the 
"boil" (i.e. at the surface, directly above the sewage 
outfall) . All other stations were located relative to 
this station. 

b. Station M-2 

Station M-2 was located on a line parallel 
to the shore through the boil and "north" along the coast 
100 to 125 meters. Both surface and bottom samples were 
taken . 

c. Station M-3 

Station M-3 was located on a line parallel 
to the shore through the boil and "south" along the coast 
100 to 125 meters . Both surface and bottom samples were 
obtained . 
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d. Station M-4 



This station was located midway between the 
shore and the boil. Both surface and bottom samples were 
taken here. 

e. Station M-5 

This station was located on a line perpendicular 
to the shore and passing through the boil, 100 to 125 m 
seaward from the boil. Both surface and bottom samples 
were taken here. 

f. Station M-6 

This station was located on the same line as 
Station M-5, but 300 m seaward from the boil. A surface 
sample only was taken here. 

g. Station M-7 

A surface sample was taken near "B" buoy 
(36-39.2' N x 121.52.9’ W) (see Pig. 2) on October 12. 

h. Stations S- and PG- 

S- stations were taken at the Seaside sewage 

outfall at the surface only. PG- stations were taken at 

the Pacific Grove sewage outfall - again only at the surface. 

i. Stations 500S to 500N 

These stations were sampled near the Monterey 
outfall in conjunction with Hopkins Marine Station. They 
lie along a line of 260.5° (true). The numbers indicate 
the number of meters the station is located either north 
or south from the boil. 
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2. On Station Procedure 



Three different vessels were utilized for sampling. 

On September 16th, October 5th, and October 12th, the NPS 
JlO-foot open motor launch was used. On September 28th, 
the NPS 63-foot research boat (a converted air-sea rescue 
vessel) was employed. The final samples were obtained in 
conjunction with Hopkins Marine Station aboard their research 
boat TAGE. 

The surface water samples for analysis were obtained by 
dipping a polyethylene bucket over the side of the vessel. 

For ease of sample withdrawal, a valve was fitted near the 
bucket's bottom. 

Sub-surface samples were taken with a portable winch 
and a messenger-tripped Van Dorn-type sampling bottle which held 
enough water for all sample analyses. The sub-surface 
samples were collected approximately one meter from the 
bottom. On November 16th, sub-surface samples were collected 
by divers from Hopkins Marine Station by manually tripping 
and returning the sampling bottle to the surface at each 
station . 

3 . Station Navigation 

On September 16th, stations were visually located in 
relation to the Monterey outfall boil. On September 28th, 
fixes were obtained by taking the bearings of known land 
marks. Fixes by this method have an error radius of 15 m. 
Beginning on October 5th, station locations were obtained 
by horizontal sextant angles (the horizontal angles between 
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three known land marks). This method of navigation yields 
a radius of error of 7 m. 

B. SAMPLE ANALYSIS 

Sample analyses were carried out immediately following 
the collection of the samples, with the exception of the 
salinity samples which may be stored for several weeks 
without deterioration. 

1 . Transmittance 
a. Equipment 

The samples for this analysis were collected in 
22-liter polyethylene bottles, approximately one-third 
filled. Analysis was performed with a modified Beckman DU-2 
spectrophotometer (Figs. 4 and 5)* The instrument was 
mounted on a 305 cm x 20.8 cm aluminum channel. The original 
cell compartment was replaced by a light-tight box, 

110 cm x 30.5 cm x 21.6 cm. Cables connecting the detector 
(D, Fig. 5) and the instrument proper (M, Fig. 5) were 
extended so that the detector could be mounted on the 
end of the light-tight box. The inside of the box was 
painted flat black and had supports to hold a glass cell, 
approximately one meter in length and an inside diameter of 
6.8 cm (C, Fig. 5). The end plates of the cell were con- 
structed of plate glass and held in place by a rubber 
0-ring and flange assembly. The total cell length from the 
inside of one cell end-plate to the other was 1.016 m. 

The cell was filled and emptied through inlet and 
outlet fittings in the cell wall. From these fittings 
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rubber hoses ran out of the light tight box to the water 
sample and vacuum aspirator. 

To reduce the effects of reflection inside the cell, 
a series of seven equally spaced baffles, having one-inch 
diameter, circular, black aperatures, was installed 
(B, Fig. 5). Control of beam divergence was accomplished 
by placing a 25-cm focal length lens in the light path 
within the instrument mounting block (L, Fig. 5). By 
decreasing the slit-width and adjusting the lens, the beam 
was kept within the horizontal limits of the photomultiplier 
tube used as a detector. The vertical beam limits were 
constrained by the cell baffle apertures and the entrance 
to the detector. This allowed variation of reflections 
in the vertical plane. When the cell was filled with water, 
the light was scattered and reflected off of the cell wall. 
With a water sample, the beam was more diffuse and covered 
a larger portion of the detector than a beam passing through 
air. Consequently, when the cell was filled with water, 
observed transmittances were relatively higher than in air. 
Absolute values were as much as 10 percent too high. In 
some cases, distilled water had transmittances of 1.10. 

Some calculations for dominant wavelength and percent 
purity were performed with "estimated correct values." When 
compared with calculations made with observed values , the 
dominant wavelengths and percent purities were in all cases 
found to agree within 5 percent. In most, the agreement 
of dominant wavelength is ±2 nm and percent purity is ±5 
percent . 
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For the reference light source , an optical light-pipe 
(P, Fig. 5) was inserted into the lamp housing and aimed 
at the tungsten filament of the lamp (T, Fig. 5). The 
other end was inserted in the side of the detector housing 
CD, Fig. 5). This end of the light-pipe was provided with 
a shutter to cut off the light and a filter to reduce the 
intensity of the reference light to within the limits of 
the photomultiplier tube. A Kodak Wratten neutral density 
filter having a density of 3.0 was used in all analyses 
except those from November 16th, when one was used having 
a density of 3*5. 

b. Transmittance Calculation 

The transmittance in air was assumed to be 1.0. 

To determine the transmittance of water at a particular 
wavelength, the ratio of the transmittance of water to air 
was taken. 

A major correction was applied to account for differences 
in reflectivity between glass-air and glass-water surfaces 
at the cell ends. The calculation was based upon the Fresnel 
reflection formula for light at normal incidence to a sur- 
face separating two media having indices of refraction n^ 
and n 2 [5, p . 179] . 

(n, - n p ) 2 

r = 

(n 1 + n 2 ) 2 

For this case the light passed through both cell ends so the 
correction had to be applied at both ends . 
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K 



where 



'wr 



wo 



ao 



K 



actual transmittance of water, 

observed transmittance of water, 

observed transmittance of air, and 

0.92 + 0.01 = a constant to account for the 
difference in indices of refraction. 



K is computed in the following manner 



K = 



1 - 2(r ) 

ga 

1 - 2 ( r J 
gw' 



where : 



and 



n - n n2 

g w 



gw 



n + n 



ga 



w 



n - n 

g a 

n^ + n o 

g a 



n g 


= 1.51 


= index 


of 


refraction 


of 


plate glass. 


n w 


= 1.33 


= index 


of 


refraction 


of 


water, and 


n a 


= 1.00 


= index 


of 


refraction 


of 


air. 


Reflection 


on the 


outside 


of 


the glass 


at both ends of the 



cell was not considered as it is the same whether air or 
water is in the cell. 

To avoid having to set the slit-width of the instrument 
such that air had a transmittance of 1.00, the transmittance 
of a water sample was obtained by taking the ratio of the 
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observed transmittance of water and observed transmittance 



of air at the same slit-width and multiplying by the 
K-factor. In this manner a "blank" of air was run for each 
day's analysis to provide reference transmittances in air 
for calculation. Transmittances for the water sample were 
observed at the same slit-widths used for the air "blank." 

Even though the reference (light pipe) side of the 
instrument was always within scale recording limits, an 
uncorrected transmittance of air or a sample of greater than 
1.0 was possible. To obtain a value within instrument range 
(0.0 to 1.05) a removable 1.0 ND filter was installed in 
the light path (F, Fig. 5) in the filter holder. This 
filter, when needed, reduced the incident light intensity 
to approximately 0.1 of the original value. The filter w as 
calibrated for each wavelength observed. To convert a 
transmittance reading with the filter to the one which 
should have been observed, the reading was divided by the 
approximate correction (approximately 0.1). 



Transmittances were observed at the following wave- 
lengths : 



375. 


.0*. 


, 400.0* 


, 410.0*, 422 


.2, 432, 


.0, 435 


.5, 438. 


.6, 444.4, 


450 . 


■ 1, 


455.9, 


461.2, 462.0, 


468.7, 


477.7, 


480. 0* : 


, 489.4, 


495. 


■2, 


510.0*, 


515.1, 529.8 


, 541.4. 


, 544.3 


, 551. 8 : 


, 561.9, 


564. 


■ 1, 


572.5, 


577.3, 584.8, 


588.7, 


599.6, 


600.8, 


610.9, 


624 , 


■ 2, 


627.3*, 


630.0*, 645. 


0, 645.5 


), 665. 


0 *, 675 . 


, 0 *, 700 . 0 *. 



These wavelengths were observed beginning on October 
12th and were not used in dominant wavelength, percent 
purity, or visual efficiency calculations. 
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2. Particle Distribution 



Water particle samples were collected in 100 ml 
polyethylene bottles. Analysis was performed on a 15-chan- 
nel Coulter counter (Model T) with a 100U diameter aperture 
The range of particle sizes which could be observed was 
between 1-04 y and 27. 6y. 

3 . Phosphate 

Samples were kept refrigerated until analyzed. The 
analysis was in accordance with the method of Strickland 
and Parsons [j 3 PP • 47-5lJ • 

4 . Salinity 

Salinity samples were collected in standard citrate 
bottles. Analysis was performed on an Industrial Instrument 
Corporation (Model RS-7B) induction salinometer. 

5 . Oxygen 

Oxygen samples were analyzed using a modified micro- 
Winkler analysis Jj4, pp . J7-J12J . 

6 . Coliform 

Analyses for coliform bacteria were performed by the 
Monterey County Health Department on the samples of October 
18th and November 16th. The results give the most probable 
number (MPN) per 100 ml sample jj3 , pp . 569 - 609 J . 

7 . Temperature 

Surface temperature was measured with a bucket ther- 
mometer to an accuracy of + 0.1°C. 



III. DATA ANALYSIS 



A. DOMINANT WAVELENGTH, PERCENT PURITY, AND VISUAL EFFICIENCY 
Corrected transmittances and chromaticity coordinates 
were calculated by a computer program. The program also 
calculated the transmissions per one-meter path-length and 
the total beam attenuation coefficients at each wavelength 
observed. Another computer program (see page 121) was 
written to graph these values as functions of their wave- 
lengths. Sample graphs are on page 117- Chromaticity 
coordinates were calculated according to the selected 
ordinate method for 30 observations [5, P- 274]. 

Tristimulus values (X, Y, and Z) were computed by 
summing the transmittances of the wavelengths listed in 
Table II and multiplying by the appropriate factor. 

TABLE II 

WAVELENGTHS (IN nm) USED FOR SELECTED ORDINATE 
CALCULATION OF TRISTIMULUS VALUES FOR STANDARD SOURCE A* 



X 


Y 


Z 


516.9 


507.7 


424.9 


561.4 


529.8 


4 36.0 


576.3 


543.7 


443.7 


587.2 


555.4 


450.5 


596.5 


566.3 


456.8 


605.2 


576.9 


462.9 


613.8 


587.9 


469-2 


623.3 


600.1 


476.8 


635.3 


615.2 


487.5 


655.9 


639.7 


508.4 


Factor: .10984 


.10000 


.03555 



Linearly interpolated from the values observed. 
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X, Y, and Z represent the amount of red, green, and blue, 
respectively, that make up the color of the sample analyzed. 
The chromaticity coefficients, x, y, and z (trichromatic 
coefficients) are the normalized tristimulus values. Per- 
cent visual efficiency is the Y tristimulus value multiplied 
by 100. With x and y, the dominant wavelength was obtained 
from a large detailed chromaticity diagram. 

•B. PARTICLE AREA AND VOLUME 

Ragweed pollen (19k - 20y in diameter) was used to deter- 
mine that channel "1.5" corresponded to a median particle 
diameter of 19 . 5y on the Coulter counter. Prom this the 
particle diameter of each of the 15 channels was found 
(see Table III) . 



TABLE III 

CHANNEL NUMBER AND PARTICLE DIAMETER (y) 



Channel 

Number 


Diameter 

(y) 


Channel 

Number 


Diameter 

(y) 


0 


27.6 


8 


4 . 34 


1 


21.9 


9 


3.45 


2 


17. 4 


10 


2.74 


3 


13.8 


11 


2.17 


4 


10.9 


12 


1.72 


5 


8.69 


13 


1.37 


6 


6.90 


14 


1.09 


7 


5.^7 







A particle area and a volume were calculated for each of 
the 15 channels on the basis that the particles were all 
perfect spheres. A computer was programmed to calculate the 
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area, volume, and totals for each of the samples. See 
page 119 for an annotated listing of this program. The 
analysis of all samples is on page 95- 

C. CORRELATION AND REGRESSION ANALYSIS 

Preliminary correlation of data was done by graphing 
pairs of variables observed. Further examination was per- 
formed with the Naval Postgraeduate School IBM 360/67 
computer. A multivariable linear regression analysis 
program (REGRE) was used. With this program any variable 
could be selected as the dependent (predicted) variable 
and any or all of the other variables could be chosen as 
the independent (predictor) variable (s). Along with F-values, 
T-values , and regression coefficients, a correlation coef- 
ficient was also calculated. The more closely this coef- 
ficient approaches ± 1.0, the greater the dependent variable 
depends upon a particular independent variable. 
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IV. DISCUSSION 



A. DATA SUMMARY AND DISCUSSION 

Appendix A is a summary of the data obtained during the 
observation program. Appendix B lists corrected transmit- 
tances and their respective attenuation coefficients as a 
function of wavelength for each sample analyzed. There are 
periodic gaps in this data due to equipment failure or 
sampling problems. The lack of complete data sets somewhat 
limited the analysis. 

Sample dominant wavelength and sample purity agree favor- 
ably with values reported by other researchers, especially 



reported by Jerlov, was based on spectral measurements 
of either downward irradiance (sensor pointed up) or 
upward irradiance (sensor pointed down). As the light from the 
energy source (the sun) becomes less attenuated (observations 
in shallower water) the dominant wavelength increases. The 
present study was based on a 1-meter path length exclusively: 
consequently the values are lower than reported in Jerlov. 

The purities were lower and the dominant wavelengths longer 
also, because the samples analyzed were from waters in the 
near-shore region. 

The percent purity of all bottom samples, except 500S 
and 175N (November l6th), is less than 20 percent. Figure 6, 
a chromaticity diagram with sample coordinates plotted. 



in the area of the Baltic Sea 




This data 
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shows a broad range of values of dominant wavelength and 
purity exist in the areas sampled. 

The salinity is 33*5 ± .5 o/oo for almost all the 

stations regardless of location. There is much mixing going 
on near the shore where the samples were taken. The mixing 
is further evidenced by the randomness of particle counts. 
Heavy particle loads (in the range of 1.09 to 27. 6y) appear 
frequently in either the surface or bottom samples (see 
Fig. 7). 

Figure 8 relates projected particle area per 2 ml sample 
and percent purity. As projected particle area increases 
the percent purity also increases. One explanation is that 
dominant wavelength is attributable to mechanisms other 
than scattering. An increase in particle area seems to 
scatter energy at all wavelengths. The dominant wavelength 
is affected less than the other wavelengths and its purity 
is enhanced. Thus as projected particle area increases so 
does percent purity (of the dominant wavelength). Jerlov, 

[ 2 , pp . l44-150~j refers to the idea that Gelbstoff is largely 
responsible for color change - especially at longer wave- 
lengths in turbid water. This is further evidenced in 
Figure 9- Dominant wavelength does not seem highly depen- 
dent upon projected particle area. Regardless of the 
station most dominant wavelength values fall between 5^0 nm 
and 585 nm. 

The value of the maximum transmittance is very dependent 
upon projected particle area. Figure 10 is a plot of 
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projected particle area and maximum transmissions. The 
curve fits the equation 

logloCA x TO" 6 ) = .5^4 - .63^(max T) 

where : 

max T = maximum transmittance (m - "'") 

2 

A = projected particle area (y ) 

The correlation coefficient for this fit is -0.59* This 
correlation seems to substantiate the observations of Burt 
[9-1 whose work indicates a major factor in maximum trans- 
mittance is projected particle area. 

The most spacially homogeneous set of data was obtained 
on November 16th. There is little variation in projected 
particle area, particle volume, dominant wavelength, 
percent purity or salinity. On this date the wind was 15 
to 20 knots and the swell was 2 - 5 ft . The day was also 
significant in that the Monterey sewage plant chlorinator 
was shut-down part of the day. Figure 11 demonstrates 
the effect of this. Early in the day coliform was low. Ey 
1030, the coliform count was high to the north. At 1250, 
when the boil was sampled, there was an extremely high 
coliform count at the outfall. Since the effluent initially 
rises to the surface above the outfall, the effluent flowed 
out some distance and by noon had worked back and collected 
at the base of the outfall. By afternoon the "pulse" of 
coliform-laden water had moved south, carried by long-shore 
currents, to yield high values of coliform around 200 m south 
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of the boil. The phosphate curve gives similar results. 
Notice, though, that the phosphate at the base of the out- 
fall is considerably lower than the surface sample. 

Coliform and phosphate seem to be proportional (Pig. 11). 
Note that most of these data were obtained with a chlorina- 
tor off. Monterey County Health records show that the 
coliform count around the Monterey outfall is generally quite 
low (less than 50 MPN per 100 ml sample). In those cases 
this relationship would not hold. 

From Figure 13 it can be seen for the Monterey outfall 
that, excluding the boil (Station M-l), the highest concen- 
tration of phosphate is invariably found at Station M-5 , 
the seaward station. Trumbauer [6] speaks of a coliform 
concentration gradient on the seaward side of the outfall. 
This has been associated with a sewage concentration 
gradient. The high phosphate at Station M-5 seems to 
reinforce this concept. 

Visual efficiency is observed to be inversely related 
to projected particle area per 2 ml sample (see Fig. 14). 

The linear correlation coefficient is -.31. Such a relation- 
ship is reasonable since the greater the particle area in 
a water sample, the greater the light scattering. The 
scattering of the light lowers the transmittance at observed 
wavelengths and the visual efficiency is lowered. 

B. FOREL-ULE SCALE 

In order to compare a sample color to the Forel-Ule scale 
in a quantitative manner, the comparison was made with 
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dominant wavelength and purity. The dominant wavelengths 
and purities of the Forel-Ule colors were not available. 

To obtain these data, solutions of blue (copper sulfate), 
yellow (potassium chromate), and brown (cobalt sulfate) 
were mixed according to Sower [3] • Forel-Ule solutions 
were made by mixing the appropriate quantities of each 
solution. Apparent color ranged from blue (scale number I) 
through green (scale number XI), to brown (scale number 
XXI). Transmittances were measured in a 1-cm cell with a 
Bausch and Lomb Spectronic 100 spectrophotometer to obtain 
transmittances at the various wavelengths. As with the 
Beckman instrument, a correction was made for the differ- 
ences in the losses due to reflection at the glass-air and 
glass-water surfaces. For the Pyrex cells used, the value 
of K was 0.93 ± 0.01. 

Chromaticity coordinates were calculated using the 
previously described programs. In addition, plots were 
made (by computer) of wavelength and transmittance. Example 
plots are shown on page 115- Appendix C lists transmittance 
as a function of wavelength for each of the 21 scale numbers. 
Figure 15 shows the correlation of Forel-Ule numbers, 
percent purity, and dominant wavelength. 

Approximately 85 percent of the samples analyzed have 
dominant wavelengths between 570 and 585 nm. This corres- 
ponds to Forel-Ule numbers ranging from XI to XVI and is 
indicative of water with an apparent brownish color. The 
samples which fall in the "bluer" region do not correlate 
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with geographic location. One (October 12th, Station M-7) 
is located at "B" buoy, the most seaward station sampled. 

The others are found at near-shore stations - either at the 
surface or bottom. 

C. MULTIVARIABLE LINEAR REGRESSION ANALYSIS 

In order to analyze the data without filling blank data 
spaces, analysis was generally performed on either pairs 
of variables or sets of observations which had the same 
empty spaces. Three data sets were examined. They are 
listed, along with the resulting correlation coefficients, in 
Appendix D. 

Linear correlation coefficients indicate dominant 
wavelength is most closely associated with temperature 
(correlation coefficient equals .64) and to lesser degrees 
to phosphate and projected particle area. Percent purity 
was greatly influenced by both phosphate and projected 
particle area. 

It should be pointed out that this was only a 
linear analysis. More sophisticated analysis (especially 
using the logarithms of most of the observed variables), 
would probably yield higher correlation coefficients and 
better curve fits. 
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V. CONCLUSIONS 



The water near the Monterey sewage outfall has a dominant 
wavelength ranging from 570 nm to 585 nm and a purity of 2 
to 40 percent. On the basis of dominant wavelength, the 
water falls between XI and XVI on a Forel-Ule scale. Linear 
regression analysis correlates dominant wavelength most 
strongly with temperature and to a lesser extent with 
phosphate and projected particle area. Percent purity is 
related to both phosphate and projected particle area. 

Maximum transmittance and projected particle area for 
this study fit the curve described by the equation: 

logloCA x 10 - ^) = .5^ - .634(max T) 

where : 

max T = maximum transmittance (m "S and 

2 

A = projected particle area (y ). 

On November 16 , 1971, the Monterey sewage plant chlor- 
inator was off for a short time. Due to this, a "pulse" 
of coliform bacteria was traced providing indications of 
effluent distribution. The utilization of this technique 
could be used in other studies. 

Dominant wavelength, percent purity, and particle area 
or volume seem to be poor signatures of sewage effluent 
from the standpoint of observing (visually) effluent dis- 
tributions and concentrations. These indicators seem 
insensitive to the addition of effluents from sewage outfalls. 
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Figure 1. Chromaticily Diagram 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 



TABLE I 



Key to Station Locations for Figures 2 and 3* 



Date 


Station 


Location 

Latitude x Longitude or direction (°true) 
and distance (m) from the boil. 


Sept 16 


M-l 


Monterey 


boil (36°-36.37 , Nx121°-52.3'W) 


Sept 16 


M-2 


053.5° 


360 m from M-l 


Sept 16 


M-3 


242.0° 


230 m from M-l 


Sept 16 


M-4 


157.6° 


180 m from M-l 


Sept 16 


M-5 


331.0° 


290 m from M-l 


Sept 28 


M-l 


Monterey 


boil ( see Sept 16th, M-l) 


Sept 28 


M-2 


057 . 0 ° 


370 m from M-l 


Sept 28 


M-3 


241.5° 


200 m from M-l 


Sept 28 


M-4 


141.5° 


190 m from M-l 


Sept 28 


M-5 


332.0° 


280 m from M-l 


Oct 5 


M-l 


Monterey 


boil (see Sept l6th, M-l) 


Oct 5 


M-2 


135.5° 


370 m from M-l 


Oct 5 


M-3 


235.5° 


220 m from M-l 


u ct 5 


M-4 


142.0° 


150 m from M-l 


Oct 5 


M-5 


333.5° 


310 m from M-l 


Oct 5 


M-6 


330.0° 


740 m from M-l 


Oct 12 


M-l 


Monterey 


boil (see Sept l6th,M-l) 


Oct 12 


M-2 


041.0° 


310 m from M-l 


Oct 12 


M-3 


257.5° 


220 m from M-l 


Oct 12 


M-4 


0 

171.5 


140 m from M-l 
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21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 



Oct 


12 


M-5 


336.0° 


250 m fro- M-l 


Oct 


12 


M -7 


"3" Buoy (36°-39,17 'Nx121 0 -52.9'E) 


Oct 


18 


M-l 


Monterey Boil (see sept l6th, M-l) 


Oct 


18 


M-2 


077.0° 


390 m from M-l 


Oct 


18 


M-3 


226° 


460 m from M-l 


Oct 


18 


S-l 


Seaside 


Boil (36°-39.7'Nxl21°-51.35 


Oct 


18 


S-2 


042.0° 


300 m from S-l 


Oct 


18 


s -3 


225.5° 


600 m from s-l 


Oct 


18 


PG-1 


Pacific 


Grove Boil (36°-38.47»Nxl2l' 


Oct 


18 


PG-2 


103.0° 


1150 m from PG-1 


Oct 


18 


FG-3 


217.0° 


1080 m from PG-1 


Nov 


16 


500N 


080. 5° 


500 m from M-l 


Nov 


16 


400N 


080.5° 


400 m from M-l 


Nov 


16 


300N 


080. 5° 


300 m from M-l 


Nov 


16 


200N 


080.5° 


200 m from M-l 


Nov 


16 


100N 


080.5° 


100 m from M-l 


Nov 


16 


BOIL 


( see Sept 16th, M-l ) 


Nov 


16 


100S 


260.5° 


100 m from M-l 


Nov 


16 


200S 


260.5° 


200 m from M-l 


Nov 


16 


300s 


260.5° 


300 m from M-l 


Nov 


16 


400S 


260.5° 


400 m from M-l 


Nov 


16 


500S 


260.5° 


500 m from M-l 
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Fipure 2. Sampling - Locations (See Table I for station key). 



MONTEREY 
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Front View of Modified 
DU-2 Spectrometer 
(cell resting on top 
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Figure 5. Modified Beckman DU-2 Spectrophotometer. 
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Figure 6. Chromaticity Diagram with Sample Coordinates Plotted. 



TABLE IV 



Key to symbols for Figures 6-8 and 12 - lb. 



Symbol 


Date 


O 


September 16, 1971 


□ 


September 28, 1971 


A 


October 5> 1971 


O 


October 12, 1971 




October 18, 1971 


+ 


November 16 , 1971 


"Flag" 


Station 


sy 


M-l 


>*✓ 


M-2 


o- 


M-3 


« 


M-4 


0 

1 


M-5 


V/ 


M-6/M-7 



O 

P 



Miscellaneous 
(open) Surface 

(solid) Bottom 

Pacific Grove 



S 



Seaside 




Figure 7. Particle volume as a function of dominant wavelength 
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Fig. 8. Projecfed particle area as a function of % purify 
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Particle Area (microns^ per 2 ml sample) *10 




Figure 9. Projected particle area as a function of 
dominant wavelength 
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Particle Area (microns^ p er 2 ml sample) x 10' 
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Figure 10. 



Projected parficle area as a {unction of 
maximum fransmission 
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Figure 11. Coliform and Phosphate as a Function of Time and 
Distance from Poll (November 16) . 



